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Abstract 

The low-energy scattering of ortho positronium (Ps) by ortho Ps has 
been studied in a full quantum mechanical coupled-channel approach. In 
the singlet channel (total spin st = 0) we find S- and P-wave resonances 
at 3.35 eV (width 0.02 eV) and 5.05 eV (width 0.04 eV), respectively and 
a binding of 0.43 eV of Ps2. The scattering length for st = is 3.95 A 
and for st = 2 is 0.83 A. The small st = 2 scattering length makes the 
spin-polarized ortho Ps atoms an almost noninteracting ideal gas which 
may undergo Bose-Einstein condensation. 

Pacs Numbers: 34.85. -Hx, 03.75.Fi, 82.30. Gg 
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Of the two types of positronium (Ps), the ortho Ps (triplet) of large life 
time 142 ns, compared to that of the para Ps (singlet) of life time 0.125 ns, has 
been of great experimental interest. The preparation of ortho Ps from positron 
interaction with atoms and molecules is well under control. One can study the 
interaction of Ps, a matter-antimatter bound state, with different atoms and 
molecules. There has been precise measurement of total cross section of ortho- 
Ps scattering by H2, N2, He, Ne, Ar, C4H10, and C5H12 H]. Such processes may 
take place on the outer surface of a star and are also of astrophysical interest. 
There have also been studies on bound states of Ps with different atoms ^ as 
well as of Ps2 molecule . 

In addition to the interaction of Ps with matter, the interaction between two 
Ps atoms is also of general interest. As the anti-atom of Ps is also a Ps, one 
can study the interaction between matter-antimatter (Ps-Ps) states, where the 
constituents (Ps) are also matter-antimatter (electron-positron) bound states. 
In the recent study of hydrogen-antihydrogen interaction, the constituents are, 
however, not matter-antimatter bound states 

Apart from the theoretical interest in testing fundamental symmetries, (charge 
conjugation, parity, and time reversal) from a study of Ps-Ps interaction, the 
Bose-Einstein condensation (BEC) of a collection of spin-polarized ortho Ps 
(Pst) seems to be viable in the laboratory |||, |j. The critical temperature Tc for 
BEC of ideal bosons of mass m and density n is given by Tc = (27r^/m)(n/2.65)^/'^ 
. Hence, a small mass of Ps at a reasonable density should facilitate BEC by 
leading to a large Tc. 

Recently, there has been experimental realization of BEC in trapped atoms 
(^H, ^Li, 23Na, 87Rb, 4He*(23Si)) at temperatures of few hundred nK |^|. 
This initiated intense theoretical and numerical activities on the subject pO|. 
In actual experiment a trap was needed to increase the density of atoms to a 
desirable level. However, with the trap emerges a new energy scale huj{'^ few 
nK) in the problem where to is the angular frequency of the trap and a new 
Tc = 0.94:hujN^/^ where TV is the number of atoms. For N ^ 10^ this leads to 
a Tc of few hundred nK for trapped atoms. 

The interaction between the bosonic atoms at low temperature is propor- 
tional to the interatomic scattering length a [|o| . In most cases of trapped atoms 
above |a| is large which implies a strong interaction and hence a large devi- 
ation from an ideal gas. This makes a theoretical investigation via a strongly 
nonlinear Schrodinger equation extremely complicated |^| . 

In the total spin st — state, the collision between two unpolarized ortho 
Ps may lead to two para Ps and and their quick eventual annihilation. The Ps2 
molecule exists in this overall singlet state where both the electron and positron 
pairs also remain in singlet states. Hence, in a collection of unpolarized ortho 
Ps the long-lived ortho Ps atoms will disastrously transfer to short-lived para 
Ps, and for an efficient BEC one should use Pst |, |. For two Pst, ST = 2 and 
both the electron and positron pairs remain in triplet states. 

The small mass of Ps should enhance BEC, and 10^ Ps atoms in a cylindrical 
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cavity in Si of volume 10^^^ cm"^ should allow BEC at tens of Kelvin temperature 
according to the ideal gas BEC formula for Tc ||, ^, in contrast to hundreds 
of nK for trapped gas |j . The advantages of a Si cavity over other metalic 
cavities for cooling Ps atoms and achieve BEC have been discussed in Ref. Q. 
If n|ap << 1, the BEC is close to one of an ideal gas |5|. It is interesting to 
calculate the relevant scattering length for Ps-Ps scattering and see if the BEC 
of Pst can be considered to be ideal or nearly ideal. 

In this Letter we present a full quantum mechanical three-Ps-state coupled- 
channel method with Ps(ls)Ps(ls,2s,2p) states for the study of scattering of 
ortho Ps by ortho Ps in both st ~ ^ (singlet) and st = 2 states. We do 
not consider the simultaneous excitation of two Ps atoms in this study. The 
St = state is interesting for studying the Ps2 molecule and the sr = 2 state 
for studying the interaction among the Pst atoms to be used in BEC. Similar 
coupled-channel approach was successfully used in Ps scattering by H He 
[||, Ar, Ne ||, Li, Na, K, Rb, Cs ||, and H2 In the Ps-H system 

we could reproduce the variationally known binding and resonance energies 
and scattering length in the singlet channel. In Ps-Li, Ps-Na, Ps-K we could 
reproduce the variationally known binding energies @ in the singlet channel 
and predict resonances in S, P, and D waves near the lowest Ps excitation 
threshold. In Ps-He, Ps-Ar, Ps-Ne, and PS-II2 we could give a good account 
of the experimental total cross section at low and medium energies. Hence it 
seems worthwhile to test this approach in Ps-Ps scattering. 

Here we calculate the Ps-Ps scattering lengths, low-energy phase shifts and 
elastic and inelastic cross sections in both st = and st = 2 channels. In the 
singlet channel (st = 0) we find resonances in S and P waves below the lowest 
Ps excitation threshold. 

The theory for the coupled-channel method of Ps scattering has appeared 
in the literature ^]. The relevant coupled Lippmann-Schwinger scattering 
integral equation in momentum space is given below 

/±,(k',k)=i3±,(k',k) 

27^2 2^ 7 ^ fc2/2-g2/2 + i0 ' ^ ' 

where f^,^ is the scattering amplitude, and B^,^ is the corresponding Born 
amplitude for = and 2 respectively, v and v' denote initial and final states 
of one of the Ps while the other being maintained in the Is state, k^ is the 
on-shell relative momentum of Ps-Ps in the channel ^. We present all equations 
in atomic unit (a.u.) where Ti — e — ciQ^m^X, where e (m) is the electronic 
charge (mass) and oq is the Bohr radius. Because of internal symmetry all 
direct Born amplitudes with the excitation of only one Ps, as considered in this 
Letter, are zero and the scattering proceeds solely through exchange interaction 
by successive exchange of a pair of electrons between two positrons. 
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The exchange Born ampUtude for Ps-Ps scattering can be easily derived as 
in Ref. [|l2j and the details will be elaborated elsewhere. The time-reversal 
symmetric exchange Born amplitude is given by 

X / xL(r')e''^ ""'/'xis(r')dr', (2) 

with D = (fc2 + fc'2)/4+[/3L + (/3^+/3^')/2], where Q = k - k' is the momentum 
transfer, x(r) the Ps wave function, and and 01, are binding energies of 
initial and final states of Ps in a.u., respectively. The exchange amplitude with 
+ sign in Eq. refers to an attractive interaction and corresponds to overall 
St — ^ state mentioned above where the electron and positron pairs are in singlet 
states. The exchange amplitude with — sign refers to a repulsive interaction and 
corresponds to the overall sr = 2 state with both the positronium pairs in the 
triplet state. This exchange Born amplitude for Ps-Ps scattering is considered 
to be a symmetrized generalization of the Ochkur-Rudge exchange amplitude 
for electron scattering [T^ ]. 

After a partial- wave projection of Eq. (|l|), the resultant one-dimensional 
scattering equations are solved by the method of matrix inversion. Forty Gauss- 
Legendre points are used in the discretization of each momentum-space integral. 

First we study the binding of Ps2 using the present coupled-channel approach 
with Ps(ls)Ps(ls,2s,2p) states. In the singlet channel, this approach leads to 
a Ps2 molecule of binding 0.43 eV in excellent agreement with the accurate 
binding of 0.435445 eV |^ obtained by Frolov and Smith, Jr. using variational 
method. 

In Fig. 1 (a) we present the Ps-Ps singlet (st = 0) phase shifts at different 
Ps-Ps center of mass energies (-Ec.m. = 13.6fc^ eV) in S, P, and D waves. It 
should be noted that i?c.m. =5.1 eV corresponds to the first excitation thresh- 
old of one of the Ps. In Fig. 1 (b) we present the same for the st = 2 channel. 
From Fig. 1 (a) we find that there are S- and P-wave resonances in the singlet 
channel below the lowest excitation threshold of the Ps where the phase shifts 
jump by TT. The position and the width of the resonances are obtained by fit- 
ting the corresponding cross sections to a Breit-Wigner form. The present S- 
and P-wave resonance energies are 3.35 eV (width 0.02 eV) and 5.05 eV (width 
0.04 eV), respectively. The appearance of these resonances in the electronic 
singlet state of Ps scattering by atoms with a single active electron seems to be 
a general phenomena and similar resonances have been found in Ps scattering 
by H and alkali- metal atoms |Q. The ability of the present approach to 
reproduce the precise binding of Ps2 as well as to predict resonances in S and P 
waves assures us of its realistic nature and the legitimacy to use it in the st = 2 
channel, where one has essentially the same set of coupled equations (0) with 
a change of sign of the exchange Born amplitudes imposed by the condition of 
antisymmetrization for identical fermions. 
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Next we report the results for scattering lengths. In the singlet channel 
ao — 7.46 a.u. — 3.95 A, and in the st — '2 channel a2 = 1.56 a.u. = 0.83 A. 
The first is in agreement with the prediction of Platzman and Mills based 
on qualitative argument: ao — 3 A, and should be considered to be refinement 
over this estimate. In the interaction of two Pst, only 02 is pertinent as the 
two electrons as well as two positrons have their spins aligned parallel (triplet 
states) in two Pst. It is this scattering length which governs the BEC of Psj. 

We exhibit in Fig. 2 the present elastic Ps{ls)Ps(ls), and inelastic Ps(ls)Ps(2s) 
and Ps(ls)Ps(2p) cross sections at different cm. energies for Ps-Ps scattering 
in the two spin states: st = (full line) and 2 (dashed line). The corresponding 
exchange Born contribution (dashed-dotted line) to the cross sections is also 
shown. At high energies the cross sections for both st = and st = 2 states 
tend to the corresponding first Born results and hence to each other as the two 
Born estimates are the same. We find from Fig. 2 that the energy-dependence 
of the elastic cross section is similar to that for Ps scattering by H and 
alkali- metal atoms [|4|. The cross section has a monotonic slow decrease with 
increasing energy. In the st = channel it also has a local minimum and max- 
imum below the lowest Ps excitation threshold manifesting the S- and P-wave 
resonances. 

Let us see the consequence of the present calculation on the proposed BEC 
of Pst |, I of 10^ Pst atoms in a Si cavity of macroscopic volume 10 cm 
leading to a density n of 10^* atoms/cm^. The interaction between these atoms 
will be governed by the scattering length 02 ^ 0.8 Aand n|a2p gives a good 
measure of the nonideal nature of the condensate. In this case ri| 02]^ ^ 10^^ << 
1, and the ideal gas BEC critical temperature formula without trap should hold 
good. The smallness of the scattering length 02 for Psf makes it an excellent 
candidate for very weakly interacting BEC. We recall that for the experimentally 
observed BEC in trapped atoms, the relevant scattering length for ^^Na is 27.5 
A, for ^^Rb is 57.7 A, for ^Li is -14.5 A, for '*He*(23Si) is 16 A§ |. The larger 
scattering lengths make these atoms strongly interacting and their condensates 
should have larger deviation from ideal-gas BEC. In the case of ^Li the large 
negative scattering length implies strong attraction between the atoms which 
leads to condensate which even experience collapse as the condensate grows in 
size 1^. 

Next we study the profile of the BEC of N atoms of Psf in a spherical Si 
cavity of radius R in the ground S state described by the following mean-field 
Gross-Pitaevskii (CP) equation in atomic units |l^ 



l(f(h(r) Na' " ' ^ 



4 dr^ 



(r) 



^(r) = fi^ir), (3) 



where /i is the chemical potential and the radial wave function ?/'('') = </'('')/'" is 
normalized according to |(/)(r)pdr — 1. As proposed ||^ we consider N = 10^ 
in a spherical cavity of volume 10^^'^ cm'^ which corresponds to i? = 5440 a.u. 
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We solve Eq. (||) variationally subject to the boundary condition ^(0) — (t>{R) = 
as in a spherical box. In the ideal gas case 02 = and the normalized solution 
of Eq. (|) is (t){r) = ./2/Rsm{TTr / R) with ^l = tt'^/{AR^) = 8.34 x 10"® a.u. In 
Fig. 3 we plot the variational solution of Eq. (^) for — 10^ and R = 5440 
a.u. and compare with this ideal-gas result. The variational result for is 
4.18 X 10~^ a.u. The increase in the value of fi compared to the ideal-gas case is 
due to the repulsion among the atoms, which contributes positively. In the ideal 
gas case the center of the spherical cavity is the region of maximum density. In 
case of Pst the region of maximum density has moved away from the central 
region due to the repulsive interaction between the Pst atoms. However, as the 
repulsion among the atoms is weak the wave function of the BEC of Psj atoms 
is qualitatively similar to the ideal gas solution. For BEC of trapped bosons 
because of the attraction of the trapping potential the center of the trap is the 
region of maximum density for both attractive and repulsive atomic interactions 
[0. 

In conclusion, we have studied the Ps-Ps scattering in a coupled-channel 
framework allowing the excitation of one of the Ps atoms to Ps(2s,2p) states. 
We report the results for scattering lengths for total spin st = (3.95 A) and 
St = 2 (0.83 A), and Ps2 binding energy (0.43 eV). We also present results 
for phase shifts in lower partial waves and elastic Ps(ls)Ps(ls) and inelastic 
Ps(ls)Ps(2s,2p) cross sections for st = and st = 2. We find resonances 
in S and P waves in the singlet state. The consideration of scattering length 
shows that Pst atoms are extremely weakly interacting and constitute excellent 
candidates for nearly ideal-gas BEC in macroscopic Si cavity without magnetic 
trap at a relatively high temperature of tens of Kelvin. By explicitly solving the 
GP equation we show that the wave function of the Pst in a spherical cavity of 
volume 10"^'^ cm'^ is similar to that of ideal gas atoms in the same cavity. 

The work is supported in part by the CNPq and FAPESP of Brazil. 
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Figure Caption 

1. S- (full line), P- (dashed line), and D-wave (dashed-dotted line) phase 
shifts of Ps-Ps scattering at different cm. energies in the (a) st = and (b) 
St = 2 states. 

2. Cross section of Ps-Ps scattering at different cm. energies in the st = 
(full line) and = 2 (dashed line) states. The corresponding Born contribu- 
tions (dashed-dotted line) are also shown. 

3. The wave function V(r-) = <?i>(r)/r of the BEC of Psj atoms (full Une) 
contrasted with the ideal gas solution (dashed line) in atomic units. 
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Figure 1 (a) 
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Figure 1 (b) 
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Figure 2 
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